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1. Factorization theorem and evolution equations for gluon TMDs 

Much effort has been devoted to properly describe the transverse momentum dependent parton 
distributions (TMDs), in order to avoid undesired features preventing them from representing phys¬ 
ical hadronic quantities - such as spurious rapidity divergences. Recently, improved definitions 
were introduced, regulating and canceling these divergencies and thus removing the bad features 
from the quark TMDs [1, 2, 3]. We extend these methods and definitions to properly define also 
the gluon TMDs in [4], 

In order to properly define all the (un)polarized gluon TMDs we derive a factorization theorem, 
separating the physics at different scales, for the transverse momentum spectrum in Higgs boson 
production through gluon fusion, and explicitly verify it at next-to-leading order (NLO). The factor¬ 
ization formula is used to derive the evolution equations for all the gluon TMDs, which are driven 
by a single universal kernel. For intermediate transverse momenta of the Higgs boson, the TMDs 
can be re-factorized and matched onto the canonical PDFs. Combined with the evolution equations 
this allows the resummation of large qj logarithms. The known perturbative ingredients allow 
us to perform the resummation of the Higgs qj -spectrum at next-to-next-to-leading-logarithmic 
(NNLL) accuracy. We supplement the resummation with modeling of the non-perturbative part of 
the TMDs and study the impact of the non-perturbative physics on the Higgs transverse momentum 
distribution. Thus we provide a general framework to deal with gluon TMDs in different processes 
and account for their perturbative and non-perturbative contributions. 

The qj distribution of the Higgs boson has received much attention, both in the context of 
perturbative QCD (see for example, [5, 6]) and soft-collinear effective theory (SCET) (see, e.g., 
[7, 8, 9, 10]). TMD gluon correlators were also considered in [11, 12]. However none of the 
previous works paid attention to the cancellation of rapidity divergences in a proper definition of 
gluon TMDPDFs. We want to emphasize, that our point of focus is on the gluon TMDs, and we do 
not aim to obtain the most precise predictions for the Higgs qj spectra. 

An interesting feature of gluon TMDs, is that there are two different distributions even in an 
unpolarized proton, one for unpolarized and one for linearly polarized gluons [13, 14]. In [15] the 
authors quantified their contribution in the context of the TMD formalism, both for Higgs boson 
and C-even scalar quarkonium {%d) and Xbo) production. We extend their efforts by implementing 
the currently known perturbative ingredients to the full extent to perform the resummation at NNLL 
accuracy, providing more accurate predictions and investigating their uncertainty. 

We make use of the effective theory point of view, in which the factorization theorem is de¬ 
rived as a stepwise matching, QCD„ /=6 —> QCD„ 5 —y SCET ?7 —y SCET Aqcd . In a first step, we 
integrate out the top-quark, leading to an effective theory for QCD with tif = 5 active flavors and 
an effective ggH vertex. Subsequently the hard modes of the gluon and quark fields, with off- 
shellness at the scale of the Higgs boson mass, are integrated out and we are left with collinear 
an soft fields described by SCET 9r . In a final step, valid only for perturbatively large transverse 
momenta A qcd « l 1t « Q, the TMDs are matched onto PDLs, with an effective description 
in SCETa qcd . In each step, we obtain a matching coefficient, correcting for the differences in the 
ultraviolet region of the theories. 

The cross section can be expressed in terms of collinear and soft matrix elements, which 
contain rapidity divergencies and thus are ill-defined, but which can be carefully combined in order 
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to cancel them and give well defined gluon TMDs, see [4] for a detailed description. There are 
different methods to deal with the rapidity divergencies by combining the soft and collinear matrix 
elements. One can use a rapidity regulator, such as the A regulator [2], the rapidity regulator 
introduced in [16] or alternatively combine the integrands in order to cancel the divergencies [1]. 
We want to emphasize that the proper definition of the TMDs does not depend on the choice of the 
regulator used [17]. 

In terms of the gluon TMDs, G g n, the cross section for Higgs qj distribution takes the form 


do 


dycf-q 


= 2 Ob O) C 2 (mj ,n)H (i mjj, H ) 


{inf 


d 2 y± e iqry± 


X GgyA( x Aiy±’ S A'^A,H)Gg/Bgv(xB,y±,SB',^B,Ll) + @{qr/m H ) ■ 


(1.1) 


C t and Hinfj.ji) = \C//(—q 2 .j.i)\ 2 are the coefficients in the two first matchings in the derivation 
of the factorization theorem. x A ji = \f^e ±y , T = {my -\-qf)/s and y is the rapidity of the produced 
Higgs boson. The Born-level cross section is (Join) = my a 2 (y )/{72n{N 2 — l)sv 2 ) given in terms 
of the Higgs boson mass mu, the strong coupling a s , the square of the proton CM energy s and 
the vacuum expectation value v. Note that the gluon TMDs depend on two scales, the factorization 
scale /i and an energy scale C - related to the separation of the two TMDs in rapidity. The evolution 
of the gluon TMDs in the renormalization scale /./ is driven by the anomalous dimensions, 


d 

d\ n /i 


InG 


\pol\ 

g/A 


{xA.b ± ,S A -,^A,P) 


- r cu S p(«^(At))ln^ - ~ Y(as(fi)) 


IHMn)) 

05, (jU) 

(1.2) 


While the evolution in the rapidity scale £ 

^ln G^{xA,b ± ,SA\^A,n) = ~Dg{b T -n), (1.3) 

is controlled by a function D g containing both non-perturbative and perturbative information. Thus 
part of the evolution equations for the TMDs has to be modeled and/or measured. Combining the 
evolution in the two scales the gluon TMDs evolve as 

G^’/a (x A ,b±,S A ', £a,/i Hf) = G^ff {x A ,bj_,S A ', ^A,i, Hi) R 8 {br', %A,i, Hi, £aj, Hf) , (1-4) 

with the evolution kernel 

R g {bT,UnHi,U,f,Hf) =ex p{^ / ^7G(o:,(A),ln^')| (^) ■ ( L5 ) 

2. Re-factorization of TMDs and r/y -resummation 

The gluon TMDs for general polarization can be decomposed into different functions for 
different proton and gluon polarizations, see e.g. [13]. For intermediate transverse momenta 
A 2 q CD « <fr « Q 2 the gluon TMDs can be matched onto their collinear analogues. There are 
three functions, ff describing unpolarized gluons in an unpolarized proton, h , " describing linearly 
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polarized gluons in an unpolarized proton and g\ L describing longitudinally polarized gluons in a 
longitudinally polarized proton, which can be perturbatively generated by the canonical (leading- 
twist) PDFs. We calculate the matching coefficients at NLO for the all three functions. 

Depending on the TMD considered, the collinear functions which describe its perturbative, 
small-/? 7 , region will be different, and also the relevant Wilson coefficients. The operator product 
expansion for the three functions are given by 

ff A {x A ,b T -, £a,aO= Y f ^C f g/j {x,bp-,^ A ,p)f j/A {x A /x;iJ.) + ^{bpA q CD ), 
j=q,q,g JxA x 

lif 8/A{2) ( XA ,b T -, Ca,M)= Y [ —Cgij(x-bT;£ Aj !d)fj/ A (x A /x-,il) + ff(b T A QC D ), 

j=q,q,g JxA x 

g 8/ i A (x Aj b T ]Q:P) = Y [ -^C 8 g/j (x,b T -,U,IJ-)gj/ A {x A /x-,ii) + ^(b T AQa)), (2.1) 
j=q,q,g JxA X 

where the matching coefficients at NLO are given in [4]. Note that the TMDs for the unpolarized 
gluons and the linearly polarized gluons are both matched onto the same PDF, but the first non-zero 
order of the matching coefficient for the linearly polarized gluons is one order higher in a s . The 
gluon helicity TMD is matched onto the helicity PDF. 

We can choose to set the resummation scale either in impact parameter space or in momentum 
space. In these proceedings we only present results with the resummation performed in impact 
parameter space, see [4] for a more detailed discussion . The resummed TMDs in impact parameter 
space can be written as 

= p{£f»(a,(P).ln|)}(|) 

x Y Cg/j(x A ,bT',£o,IM))®fj/ A (x A ',fA)) , (2.2) 

j=q,q,g 

where Co ~ jtl| ar| d /./o ~ pi, = 2e~ lE /bp- The superscript Pert signifies that it is only the perturba¬ 
tive part of the TMDs valid at small bj << I / Aqcd- 

3. Higgs qr spectrum 

For large bj we need to supplement the perturbative expression with a model with parameters 
that can be extracted from experimental data. We implement a smooth cutoff that freezes the 
perturbative contribution towards large bp. 


F g/A (x A M\U,P) = F g p ^(x A ,b T -,U,P)F NP (x A ,b T -,U), ( 3 - 1 ) 

with the cutoff prescription bp {bp) = b c \J(1 —and b, determining the separation between 
the perturbative and non-perturbative regions. We parametrize the non-perturbative piece for the 

1 We refer to the perturbative or non-perturbative nature of the transverse momentum (or impact parameter) depen¬ 
dence, leaving aside the non-perturbative PDFs. 
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Figure 1: Ratio 2% for different values of the non-perturbative parameters f f(h) and Xq, at the scales for 
Higgs boson and T)b production. The curves are calculated at NNLL accuracy and for s/s = 8 TeV. 


two TMDs contributing to Higgs production as 


pfjY = cx V [-b\(X f + X Q HQ 2 /Ql))\ , F$ p = exp[-4(A ft + A e ln(e 2 /eo))] , (3-2) 

where (2o = 1 GeV. Xq is the same for both functions, since the evolution is universal among 
all (un)polarized TMDs, and therefore, their scale-dependence is the same. The non-perturbative 
model F np should be 1 for bp = 0 not to modify the perturbative result and plays an increasingly 
important role as we increase bp. The contribution of unpolarized and/or linearly polarized gluon 
distributions in unpolarized hadron-hadron collisions depends on the process under study and has 
been discussed in several works, see for example [18, 19, 20, 21, 15]. In the cases we consider, 
Higgs boson as well as rjb production, both unpolarized and linearly polarized distributions play 
a role, and we investigate their relative contribution to the cross section. We use our results to 
quantify the contribution of linearly polarized gluons, considering the following ratio 


3g(x A ,x B ,q T ]Q) = 


f d 2 bp e-brbrtifim fa ,br\Q\Q) hf g/B(2) (. x B , bp\Q 2 , Q ) 
/ d 2 bp e-^-bTff A fa ,b T \Q 2 ,Q) ff (x B ,b T ; Q 2 , Q) 


(3.3) 


which determine the relative contribution from linearly polarized and unpolarized gluons to the 
cross section, for both Higgs boson and C-even pseudoscalar bottonium production. 

Fig. 1 shows our results for the ratio 3$ at the scales for the transverse momentum distributions 
of Higgs boson and r/^,, all at NNLL accuracy. The bands are obtained by independently varying the 
scales Co and /Iq around their default value by a factor of 2, and plotting the maximum uncertainty 
for each point in qp. In order to estimate the impact on the ratio of the different non-perturbative 
parameters, we have chosen several values in a sensible range and selected some combinations 
in limiting cases. The outcome of the numerical study is clear: the lower the scale the more 
contribution we have from linearly polarized gluons, although this contribution depends on the 
value of the non-perturbative parameters. At the Higgs boson scale the effect of linearly polarized 
gluons is small, around 1-9%, making it harder to extract their non-perturbative parameters from 
experimental data. At lower scales, as in the production of v\y. their role is enhanced, from 10% up 
to 70%, and thus experimental data can better determine them. It seems plausible that their non- 
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Figure 2: Cross section at y/s = 8 TeV (left and middle) and 13 TeV for different values of the non- 
perturbative parameters with po ~ Pb on NLL and NNLL. 


perturbative parameters can be fixed in the near future by properly combining experimental data 
for different experiments and at different scales. The framework introduced in this paper, with the 
proper definition of gluon TMDs and their QCD evolution, will be crucial in order to consistently 
address different processes in terms of the same hadronic quantities and properly extract their non- 
perturbative parameters. 

After analyzing the contribution of linearly polarized gluons for rfr, and Higgs boson produc¬ 
tion in unpolarized hadron-hadron collisions, we turn to the Higgs boson transverse momentum 
distribution at the LHC. The cross section for this process, given in Eq. (1.1), for unpolarized 
protons is 


cla 

dyd 2 q i 


= 2 Go (p ) C 2 {in 2 , p )H(m 2 H , p ) 


(2tt) 


d 2 y L e iq '- y ' 


1 

X 2 L- 


fi M {xa , b T ; Ca , A*) /f /B {*b , b T ; Cs, A*) 

+ hf g/M2) (x A , b T ; £a, jU) hf g/BU) {x B ,b T -,^B,^) + @{qT/m H ) ■ 


(3.4) 


The evolution kernel suppresses the TMDs at large bp, and this effect increases with the hard 
scale Q [22, 23]. Therefore, the larger Q the less sensitive the resummed expression will be to the 
non-perturbative contributions. The first two panels of Figure 2 show the Higgs boson transverse 
momentum distribution at v /v = 8 TeV, at both NLL and NNLL for different non-perturbative 
parameters. A significant part of the uncertainties originates in the variation of the rapidity scale, 
and the bands at NLL overlap with the NNLL bands. If we compare the two panels, we see that the 
impact of the non-perturbative contribution is rather small and induce changes to the distributions of 
similar size as the uncertainty on the prediction at NNLL. The Higgs boson transverse momentum 
distribution is thus not very sensitive to those parameters. The same conclusion was drawn in [9], 
where a Gaussian model was used to parametrize the non-perturbative contributions. The rightmost 
panel of Figure 2 shows the predictions for the distribution at yfs = 13 TeV, for different values 
of the non-perturbative parameters. The cross section is larger than at = 8 TeV, but the same 
conclusions hold regarding the sensitivity to the non-perturbative parameters. It therefore seems 
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unlikely that experimental measurements of the Higgs cjj distribution at the LHC will be precise 
enough to fix the non-perturbative parameters of gluon TMDs, apart from excluding the most vivid 
parameter values. 
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